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SUMMARY 

Zsako  has  modified  a curve  fitting  method  reported  by  Doyle  and 
in  this  report  the  method  is  further  simplified  to  allow  the  use  of  a 
digital  computer.  The  best  description  of  the  thermogram  step  was  found 
to  be  given  by  19  equally  spaced  co-ordinates  between  the  (1  - a)  limits 
of  0.05  to  0.95. 

These  co-ordinates,  together  with  the  heating  rate,  are  used 
to  calculate  E and  log  Z,  by  an  iterative  process  which  minimises  the 
standard  deviation  of  the  log  g(a)  - log  p(x)  values. 

The  standard  deviation  allows  a comparative  assessment  of  the 
validity  of  the  kinetic  equations  chosen  to  represent  the  mechanism  of  the 
reaction.  A further  simplification  involved  the  use  of  an  approximate 
equation  for  - log  p(x)  in  place  of  tabulated  values. 

The  TGA2  method  of  calculation  was  shown  to  be  valid  by 
re-calculating  E and  log  Z for  standard  curves  which  were  generated  from 
tables  of  - log  p(x). 

Published  data  by  Doyle  and  Zsako  were  re-calculated  and  the 
answers  were  in  good  agreement  with  the  original  results. 


INTRODUCTION 


The  thermal  decomposition  of  solids  may  be  described  by  the 
equation  A(g)  - B(a)  + C(g) . 

This  is  a simplified  statement  of  a complex  process  in  which  many 
different  steps  may  be  proceeding  simultaneously.  These  include  physical 
processes  such  as  diffusion  of  heat  and  diffusion  of  gaseous  products,  as 
well  as  breaking  of  chemical  bonds. 

Thus,  the  kinetics  of  thermal  decomposition  reactions  depend  upon 
experimental  conditions  as  well  as  the  chemical  nature  of  the  solid  [1,2]. 

Many  efforts  have  been  made  to  obtain  kinetic  parameters  from 
isothermal  and  dynamic  thermograms  and  these  have  been  reviewed  [3,4]. 

Utilising  a single  weight  loss/temperature  curve  has  an  obvious 
advantage  over  Isothermal  weight  loss/time  curves  produced  at  several 
different  temperatures.  Graphical  comparison  methods  have  been  used  [5] 
but  if  TGA  data  can  be  prepared  in  a suitable  form,  the  computer  can 
remove  most  of  the  tedium  and  cost  of  the  calculation. 


A suitable  method  for  this  kind  of  treatment  was  a curve 
fitting  method  reported  by  Doyle  [6]  and  later  simplified  by  Zsako  [7] 
This  report  describes  how  the  method  was  further  simplified  to  allow 
the  required  parameters  to  be  extracted  quickly  while  retaining  the 
accuracy  of  the  original  method. 
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The  usual  notations  were  used,  ies- 
a “ fractional  weight  loss, 

E - activation  energy  in  kcal/mole, 

Z ■ pre-exponential  factor, 
b " order  of  reaction, 

R - gas  constant, 

S*  ■ activation  entropy, 
h “ Planck '8  constant, 
k » Boltzmann’s  constant, 

Tj  ■ temperature  of  half  weight  loss, 
u - E/RT, 

x - the  value  of  u at  the  sample  temperature. 
2.  THE  DOYLE/ZSAKO  METHOD 

2.1  The  outline  of  the  method 


The  equation  to  a thermogram  is 


g(o)  ■ -^POO*  ^ 

where  the  form  of  g(a)  depends  upon  the  nature  of  the  process  and 
-x  00  -u 

p(x)  ■ - I - — .du.  (2) 

x _ u 


Doyle  prepared  a table  of  - log  p(x)  from  x “ 10  to  x ■ 50.  A 
first  approximation  of  E was  calculated  assuming  the  reaction  order  to  be 
unity,  and  this  approximation  was  improved  by  trial  and  error  curve 
fitting  techniques.  Zsako  simplified  this  process  as  follows.  Taking  logs 
and  rearranging,  equation  (1)  becomes 

log  ■ log  g(a)  - log  p(x>  - B,  .....(3) 

where  B is  dependent  upon  the  nature  of  the  compound  and  the  heating  rate 
but  not  upon  the  temperature. 

The  value  of  g(a)  can  be  found  for  a given  temperature  if  the 
particular  function  of  a is  known,  and  - log  p(x)  can  be  obtained  from 
Doyle's  tables  if  E is  known. 


...  ... 


A 


<a 


Then,  if  E ia  known  and  used  to  calculate  B at  several  temperatures, 
the  constancy  of  B will  indicate  the  validity  of  the  particular  kinetic 
equation  selected  to  calculate  g(a).  Conversely,  if  the  kinetic  equation 
representing  the  reaction  is  known,  the  constancy  of  B will  indicate  the 
validity  of  a chosen  E value.  The  constancy  of  B can  be  quantified  by 
means  of  the  standard  deviation. 


SIGMA 


A*  - 

/ N - 1 


(4) 


where  N is  the  number  of  B values.  (SIGMA)  will  indicate  which  of  two 
E values  (or  which  of  two  kinetic  equatlonsj^est  fits  the  reaction. 

Taking  the  form  of  the  kinetic  equation  as  f(a)  ■ (1  - o)^,  Zsako 
assumed,  in  turn,  reaction  orders  of  0,  1/3,  1/2,  2/3,  1 and  2 and  found 
a value  for  E,  in  each  case,  which  gave  a minimum  standard  deviation  of 
B values.  Comparing  the  (SIGMA)  values  corresponding  to  the  six  E values, 

the  overall  (SIGMA)  value  indicated  which  of  the  six  kinetic  equations 
best  fitted  the  reaction. 


(5) 


f] 

f 1 


f 
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Then  1 ■ log  ZE/Rq  and 
log  Z • B + log  Rq  - log  E, 
also  tha  apparent  activation  entropy 

S*  ■ 2.303R  log  rjr  • (6) 

i 

In  order  to  extend  Doyle's  tabulated  - log  p(x)  values  when 
x > 50  Zaako  used  the  approximation 

P(x)  - e‘x(^-  - (7) 

and  Intermediate  values  of  - log  p(x)  were  found  by  interpolation. 

2.2  Modification  of  the  Doyle/Zsako  method  for  TGA2 

2.2.1  Calculation  of  - log  p(x) 

Comparing  the  approximation  equation  (7)  with  values  for  x 
obtained  from  standard  tables  of  e~x  and  the  Integral 


“ -u 


it  was  shown  (table  2)  that  the  approximation  is  sufficiently  close  down 
to  x • 4.0.  Then,  instead  of  using  tables  of  - log  p(x),  equation  (7) 
was  used  throughout  this  work.  In  practice  x is  rarely  less  than  4.0, 
but  in  the  event,  a diagnostic  message  was  printed  by  the  computer 
program  to  Indicate  that  E was  inaccurate. 

TABLE  2 

A Test  of  the  Approximation  Equation  (7) 
over  the  Range  x - 2 to  50 


X 

A 

- Log  p(x)  from  Integral 
(Tabulated  Values) [8] 

B 

- Log  p(x)  from 
App  roximat ion 
Equation  (7) 

Ratio  B/A 

50 

25.129 

25.125865 

0.99988 

40 

20.597 

20.594464 

0.99988 

30 

16.010 

16.010155 

1.0000 

20 

11.328 

11.331665 

1.0003 

10 

6.415666 

6.438694 

1.00359 

9 

5.896899 

5.925212 

1.00480 

8 

5.368888 

5. 404 500 

1.00663 

7 

4.829246 

4.875506 

1.00958 

6 

4.274602 

4.337379 

1.01469 

5 

3.699839 

3.790578 

1.02453 

4 

3.096592 

3.241743 

1.04687 

3 

2.449659 

2.733754 

1.11597 

2 

1.726291 

16.882666 

9.780 
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3.  DESCRIPTION  OF  THE  PROGRAM  TGA2 

TGA2  is  written  in  FORTRAN  and  has  been  successfully  used  for 
several  years,  and  recently  on  the  IBM  370/168.  It  is  described  under  the 
headings  Input,  Program  and  Output,  in  sections  3.1  to  3.3.  A thermogram 
may  consist  of  one  or  more  steps  which  represent  different  reactions  and 
these  steps  must  be  treated  separately.  The  data  for  a weight  loss  step 
were  plotted  as  active  weight  fractions  (1  - a),  as  a function  of  temperature 
and  a smooth  curve  was  drawn  using  a flex  rule.  It  was  the  position  on  the 
temperature  scale  and  the  shape  of  this  smoothed  curve  which  was  used  as 
a basis  for  the  calculation  of  kinetic  parameters.  Values  of  temperature 
were  taken  from  the  smoothed  curve  corresponding  to  equal  increments  of 
(1  - a)  and  an  example  is  given  in  figure  1 where  19  co-ordinates  were 
taken  corresponding  to  (1  — a)  values  from  0.05  to  0.95  in  equal 
increments  of  0.05.  The  extremes  of  the  curve,  ie,  (1  - a)  - 1.0  and  0.0 
were  not  used  because  of  the  inaccuracy  of  estimating  temperature  at  small 
slopes. 


3.1  Input 


Line 

Number 

Format 

Symbol 

Explanation 

1 

12 

NDATA 

The  number  of  pairs  of  co-ordinates 
describing  the  smoothed  curve 

2 

F6.3 

HDIF 

Increments  of  (1  - a) 

3 

6A4 

ORDER(J) 

Titles  for  the  different  kinetic 
equations 

4 to 

(NDATA  + 4) 

F7.4 

H(I) 

Active  weight  fractions  (1  - u) 

NDATA  + 5 

18 

RUN 

Run  identity  number 

A8 

STEP 

Step  identity  number 

A8 

RANGE 

Temperature  limits  of  step  in  *C 

F8.3 

QM 

Heating  rate  ”C/min 

NDATA  + 6 to 

F7.1 

TC(I) 

Temperature  in  "C  corresponding 

(2  NDATA  + 6) 

to  the  active  weight  fractions  on 
line  4 onwards 

Last  line 

_ 

— 

Blank 

Note  that  several  sets  ot  data  may  be  calculated  in  on*r. 
computer  run,  additional  sets  corresponding  to  lines  (NDATA  + 5)  to 
(2NDAIA  + 6)  being  stacked  together  without  spaces.  In  all  c*ses  the 
last  line  of  the  complete  data  input  must  be  blank.  All  data  sets  applied 
at  one  time  must  be  compatible  with  the  lines  1 to  (NDATA  + A)  which  are 
not  repeated. 


3.2 


3.2  Program  TGA2 

A flow  chart  showing  the  logic  of  TGA2  la  shown  In  figure  5 and 
a full  prlnt-out  of  the  FORTRAN  statements  Is  reproduced  In  appendix  A. 
The  main  calculations  embodied  In  TGA2  are  explained  In  sections  3.2.1 
to  3.2.3. 


3.2.1 


Minimum  of  SIGMA  (equation  (4)) 


At  a constant  heating  rate,  SIGMA  is  a function  dependent  on 
activation  energy.  Then,  for  a given  rate  equation,  g(a) , SIGMA  may  be 
minimised  by  finding  the  appropriate  value  for  E.  Subroutine  VD01A  [9] 
finds  the  minimum  of  a function  of  a single  variable  by  an  iterative 
process.  In  this  case  the  variable  is  E and  the  function  is 

/5H? 

/ N - 1 

The  subroutine  requires  a starting  approximation  of  E and  this 
is  provided  by 

E,  . • — — — 7.  — ^ ~ (cf,  reference  [10], 

(approx)  - *n(l  - a)  dT 

equation  21)  .....(8) 


equation  21)  .....(8) 

The  iterative  process  continued  until  E is  precise  to  within 
25  cal /mole. 


3.2.2  Calculation  of  — — L 

Subroutine  SLOPE  calculates  the  slope  of  a curve  from  a Hat  of 
co-ordinates  using  the  Newton-Stir ling  [11]  formula  for  central  differences 
when  the  argument  changes  by  a fixed  amount  (HDIF  in  this  subroutine) . 

The  slope  at  (!-<*)■  0.8  is  used  in  equation  (8)  to  calculate 


3.2.3 


Diagnostics 


If  E/RT  is  allowed  to  assume  a value  < 2.0,  the  program  will 
stop  at  Statement  117  having  tried  to  find  the  log  of  a negative  number. 
To  avoid  this  during  the  iterative  process,  a barrier  has  been  inserted 
(Statement  109)  so  that,  if  E/RT  < 2.1,  SIGMA  is  put  to  1000.0  to  change 
the  direction  of  search  for  SIGMA^^. 

If  E/RT  is  < 4.0  for  SIGMA , , table  2 shows  that  E will  be 
inaccurate  and  in  the  event  a statemSn?  is  printed  to  that  effect. 

If  the  number  of  evaluations  of  SIGMA  exceeds  the  limit  stated 
when  VD01A  is  called  (100  in  TGA2)  then  a diagnostic  is  printed  - 
"MAXFUN  TOO  SMALL". 

A diagnostic  message  will  be  printed  if  ROUNDING  ERRORS  prevent 
the  calculation  of  E to  the  desired  precision. 


-u 


3.3 


Output 


The  output  from  TGA2  ia  in  the  form  of  a table  and  an  example 

la  given  in  figure  2.  This  table  liata  SIGMA  and  the  correaponding 

calculated  values  for  activation  energy,  logz?  Z and  activation  entropy 

(S*) . Each  line  represents  a different  g(a)  function  and  the  best  values 

are  those  in  line  with  the  minimum  SIGMA  . (sic)  value. 

min 

Points  to  note  are:- 

(a)  The  calculation  will  indicate  the  best  fit  to  the  data 
of  the  g(a)  functions  included  in  the  program.  If  the  correct 
g(a)  function  is  not  included,  a misleading  answer  will  be 
given. 

(b)  Log  Z and  S are  calculated  in  all  cases,  but  these  are 
valid  only  if  SIGMAm<n  is  sufficiently  small. 

4.  DESCRIPTION  OF  THE  CURVE 

The  basis  for  the  calculation  of  kinetic  parameters  by  TGA2 
was  a description  of  the  smoothed  thermogram  step  in  terms  of  the 
co-ofdlnates  (1  - a)  and  T.  There  were  possible  variations  in  this 
description,  such  as  range  and  number  of  co-ordinates,  which  would  possibly 
affect  the  results  of  the  calculation.  These  variati-  is  were  investigated 
as  follows 

(a)  The  effect  of  limiting  the  range  of  (1  - o)  by  truncating 
both  ends  of  the  curve  simultaneously  (section  4.2). 

(b)  The  effect  of  limiting  the  number  of  co-ordinates, 
uniformly  spaced  with  respect  to  (1  - a)  (section  4.3). 

(c)  The  effect  of  truncating  the  ends  of  the  curve  separately 
(section  4.4). 

(d)  With  a (1  - o)  spacing  of  0.05  (19  co-ordinates)  the 
effect  of  taking  the  top  7,  middle  7 and  bottom  7 co-ordinates 
in  turn  (section  4.5). 

In  order  to  make  a quantitative  assessment  of  the  various 
treatments  discussed  in  sections  4.2  to  4.5,  standard  thermograms  were 
required  with  known  kinetic  parameters.  These  were  calculated  from  assumed 
values  of  the  relevant  kinetic  parameters  and  the  equation  to  the 
thermogram  (equation  (1)).  They  covered  the  range  of  activation  energies 
which  were  likely  to  be  encountered  in  thermal  degradation  reactions,  ie, 
from  10  to  66  kcal/mole  (section  4.1). 
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It  can  be  seen  from  the  results  of  the  Investigations  into 
the  description  of  the  curves  in  cables  5,  6 and  7 that  the  values  of  E 
and  log  Z are  not  greatly  affected  by  the  variations  applied  and, 
therefore,  a different  quantitative  assessment  was  devised  depending 
upon  Che  SIGMA  values.  It  may  be  said  that  the  choice  of  best  SIGMA 
is  easier  and  probably  more  reliable  if  it  is  much  smaller  than  the  next 
best  value  of  SIGMA,  or,  in  other  words,  if  the  ratio  (next  best  SIGMA/best 
SIGMA)  is  large.  The  value  of  this  ratio  has  been  used  thoughout  this 
work  as  a criterion  when  judging  between  the  different  treatments  applied 
to  the  standard  curves. 

4.1  Preparation  of  four  standard  curves 

TABLE  3 


Assumed  Kinetic  Parameters  for  the  Standard  Thermograms 


E. 

kcal/mole 

Log  Z 

Reaction 
Order  (b) 

Heating  Rate, 
q,  °C/min 

10.0 

2.4 

1 

10 

30.0 

10.0 

1 

10 

50.0 

15.0 

1 

10 

66.0 

19.0 

1 

10 

The  values  shown  in  table  3 were  substituted  in  a rearranged 
equation  (1)  when  g(a)  ■ (1  - a)b,  ie, 

/i  \b  ZE  , . 

(1  - a) 

and  using  Zsako's  tables  for  - log  p(x)  at  different  temperatures, 
several  (1  - a)  and  T co-ordinates  were  calculated  and  used  as  a basis 
for  four  standard  curves.  The  curves  were  smoothed  and  Initially  described 
by  33  co-ordinates  with  a range  of  (1  - a)  from  0.02  to  0.98  in  equal 
steps  of  0.03.  These  are  reproduced  in  table  4. 
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TABLE  4 


Co-ordinates  Taken  from  Four  Smoothed  Standard  Curves 


(a) 

E - 10.0  log  Z - 2.4 

(b) 

E - 30.0  log  Z - 10.0 

(c) 

E - 50.0  log  Z - 15.0 

(d) 

E - 66.0  log  Z - 19.0 

(1  " o) 

Temperature,  ®C  1 

— 

(a)  10/102*4 

— 

(b)  30/10 10 

(c)  50/ 10 15 

(d)  66/1019 

0.98 

80.0* 

208.5 

314.2 

367.0 

0.95 

100.0 

221.6 

327.3 

377.9 

0.92 

112.3 

229.3 

334.4 

384.0 

0.89 

121.3 

234.5 

339.1 

388.2 

0.86 

128.5 

238.8 

342.6 

391.4 

0.83 

134.3 

242.5 

345.6 

394.0 

0.80 

139.4 

245.5 

348.2 

396.5 

0.77 

144.3 

248.2 

398.6 

0.74 

148.6 

250.5 

352.5 

mwm 

0.71 

152.7 

252.4 

354.4 

0.68 

156.3 

254.3 

356.2 

0.65 

159.6 

256.1 

357.8 

405.1 

0.62 

162.9 

257.9 

359.5 

0.59 

166.2 

259.7 

0.56 

169.4 

261.5 

362.1 

iil 

0.53 

172.5 

263.2 

363.6 

410.5 

0.50 

175.5 

264.7 

365.1 

411.7 

0.47 

178.6 

266.2 

412.9 

0.44 

181.5 

267.8 

367.8 

414.0 

0.41 

184.4 

269.3 

369.0 

415.0 

0.38 

187.4 

270.9 

370.3 

416.2 

0.35 

190.3 

272.4 

371.5 

417.3 

0.32 

193.2 

273.9 

372.6 

418.5 

0.29 

196.2 

275.3 

373.8 

419.7 

0.26 

199.4 

276.8 

375.1 

420.9 

0.23 

202.7 

278.5 

376.5 

422.3 

| 0.20 

206.2 

280.2 

378.0 

423.5 

i 0.17 

210.0 

282.0 

424.9 

) 0.14 

214.1 

284.1 

381.5 

426.4 

: 0.11 

219.0 

286.5 

383.5 

427.8 

0.08 

224.5 

289.2 

385.7 

429.5 

0.05 

231.4* 

292.5 

388.3 

431.6* 

0.02 

245.2* 

298.0 

392.8 

435.3* 

r 

^Calculated  from  p(x)  ■ e X(~2'  " “3")  • 


The  effect  of  Halting  the  range  of  (1  - a)  by  truncatin 
Che  ends  of  Che  curve 


Table  5 shows  chaC  Che  re-calculaced  values  for  E and  log  Z 
are  alcered  11 tele  by  Chls  treatment  of  Che  standard  curves,  and  they 
are  ac  worse  within  a few  per  cent  of  chose  assumed.  The  lowest  range, 
0.2/0. 8,  shows  a low  value  for  Che  ratio  (next  best/best  SIGMA)  in  all 
four  cases  and  chls  shows  chat  Che  ends  of  Che  curve  beyond  0.2/0. 8 are 
necessary  if  Che  correct  or  minimum  SIGMA  value  is  to  be  distinguished 
easily.  On  balance  the  range  0.05/0.95  is  to  be  preferred  with 
two  Instances  of  Che  highest  ratio  favouring  this  range. 

TABLE  5 

The  Effect  of  Limiting  the  Range  of  (1  - o)  Values 


sing  Co-ordinates  from  the  Standard  Thermograms 


E/Log  Z 

Range 

NDATA* 

E 

Log  Z 

Best 

SIGMA 

Next 

Best 

SIGMA 

Ratio 

Next  Best 
SIGMA/Best 

10/2.4 

0.02/0.98 

33 

10.1 

2.5 

0.00580 

0.0355 

6.12 

0.05/0.95 

31 

10.0 

2.4 

0.00187 

0.0224 

11.98 

0.08/0.92 

29 

10.0 

2.4 

0.00167 

0.0176 

10.54 

0.20/0.80 

21 

10.1 

2.4 

0.00152 

0.00693 

4.56 

30/10.0 

0.02/0.98 

33 

30.1 

10.0 

0.00525 

0.0353 

6.72 

0.05/0.95 

31 

30.0 

10.0 

0.00452 

0.0244 

5.4 

0.08/0.92 

29 

30.0 

10.0 

0.00446 

0.0194 

4.35 

0.20/0.80 

21 

30.1 

10.1 

0.00429 

0.00947 

2.21 

50/15.0 

0.02/0.98 

33 

50.2 

15.1 

0.00624 

0.0320 

5.13 

0.05/0.95 

31 

50.4 

15.1 

0.00468 

0.0250 

5.34 

0.08/0.92 

29 

50.6 

15.2 

0.00419 

0.0198 

4.73 

0.20/0.80 

21 

50.7 

15.2 

0.00355 

0.00624 

1.76 

66/19.0 

0.02/0.98 

33 

66.0 

19.0 

0.00633 

0.0312 

4.93 

0.05/0.95 

31 

65.7 

18.9 

0.00486 

0.0222 

4.57 

0.08/0.92 

29 

65.5 

18.8 

0.00356 

0.0182 

5.11 

0.20/0.80 

21 

65.2 

18.7 

0.00248 

0.00829 

3.34 

1 


Note:  In  all  cases  the  best  fit  was  given  with  b ■ 1. 
^Number  of  co-ordinates. 


4.3  The  effect  of  limiting  the  number  of  uniformly  spaced  (1  - a 


values  in  the  range  0.05  to  0.95 


Again  the  recalculated  E and  log  Z values  vary  little  with  the 
different  treatments. 
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Comparing  the  recalculated  values  of  E and  log  Z in  table  6 
with  the  assumed  values,  the  smallest  NOATA  is  marginally  favoured.  The 
ratio  (next  best/best  SIGMA)  suggests  that  the  best  distinction  is 
obtained  at  either  NDATA  ■ 7 or  19. 


TABLE  6 


The  Effect  of  Limiting  the  Number*  of  Uniformly  Spaced  (1  - q)  Values 


E/Log  Z 

(1  - a) 
Limits 

1 

(NDATA)* 

(1  - a) 
Spacing 
(HDIF) 

— 

E 

Log  Z 

1 

Next 

Best  I „ _ 

SlCHA  isiGMA 

I 

1 

Ratio 

Next  Best 

SIGMA/ Best 

10/2.4 

0.05/0.95 

31 

0.03 

10.0 

2.4 

j 

0.00187  10.0224 

12.0 

19 

0.05 

10.0 

2.4 

0.00219 

0.0243 

11.1 

7 

0.15 

10.0 

2.4 

0.00165 

0.0329 

19.9 

30/10.0 

0.05/0.95 

31 

0.03 

30.0 

10.0 

0.00452 

0.0244 

5.4 

19 

0.05 

29.9 

10.0 

0.00375 

0.0259 

6.91 

7 

0.15 

29.9 

10.0 

0.00521 

0.0353 

6.78 

50/15.0 

0.05/0.95 

31 

0.03 

50.4 

15.1 

0.00468 

0.0250 

5.3 

19 

0.05 

50.4 

15.1 

0.00278 

0.0254 

9.1 

7 

0.15 

50.3 

15.1 

0.00309 

0.0335 

10.8 

66/19.0 

0.05/0.95 

31 

0.03 

65.7 

18.9 

0.00486 

0.0222 

4.4 

19 

0.05 

65.8 

18.9 

0.00520 

0.0240 

4.6 

7 

0.15 

66.1 

19.0 

0.00676 

0.0306 

4.5 

*Number  of  co-ordinates. 


4.4  The  effect  of  progressive  removal  of  the  ends  of  the  curve 

Figures  3 and  4 show  a similar  effect  and  there  is  a downward 
trend  in  the  ratio  (next  best/best  SIGMA)  as  the  end  of  the  curve  is 
reduced.  The  maximum  at  0.05  and  0.95  in  some  cases  suggests  that  the 
co-ordinates  at  the  extremes  of  the  curve,  ie,  0.02  and  0.98,  are  not 
always  reliable.  Limits  of  0.05  and  0.95  are  favoured  on  balance. 

4.5  The  relative  importance  of  the  different  parts  of  the  (1  - q)/T 

curve 

Taking  a curve  represented  by  19  co-ordinates  with  a (1  - a) 
spacing  of  0.05,  the  re-calculated  results  from  the  whole  curve  is 
compared  in  table  7 with  the  re-calculated  values  from  the  top,  middle 
and  bottom  parts  ot  the  curve  respectively.  Surprisingly,  the  E and 
log  Z values  do  not  vary  much,  except  in  the  case  of  the  bottom  part 
of  the  66/19.0  curve.  The  list  of  ratios  (next  best/best  SIGMA)  shows 
that,  as  far  as  distinguishing  one  SIGMA  from  another,  the  middle  and  ' ne 
top  parts  of  the  curve  vere  poor,  with  an  improvement  shown  from  the 
bottom  part  of  the  curve  which  was,  therefore,  judged  to  be  the  most 
Important  part  for  this  purpose. 
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TABLE  7 


The  Relative  Importance  of  the  Different  Parts  of  the  (1  - a)/T  Curve 
with  all  (1  - a)  Spacings  of  0.05 


E/Log  Z 

(1  - a) 
Limits 

Part  of 
Curve 

Order 

E 

Log  Z 

Best 

SIGMA 

Next 

Best 

SIGMA 

Ratio 

Next  Best 
SIGMA/Best 

10/2.4 

0.05/0.95 

Whole 

i 

10.0 

2.4 

0.00219 

0.0243 

11.1 

0.05/0.35 

Bottom 

i 

10.1 

2.5 

0.00138 

0.00594 

4.3 

0.35/0.65 

Middle 

i 

10.0 

2.4 

0.000475 

0.00193 

4.1 

0.65/0.95 

Top 

2/3 

9.7 

2.2 

0.00253 

0.00261 

1.0 

30/10.0 

0.05/0.95 

Whole 

1 

29.9 

10.0 

0.00375 

0.0259 

6.9 

Bottom 

1 

29.8 

9.9 

0.00197 

0.00726 

3.7 

Middle 

1 

29.9 

9.7 

0.00115 

0.00142 

1.2 

Top 

J 

28.7 

9.4 

0.00188 

0.00208 

1.1 

50/15.0 

0.05/0.95 

Whole 

1 

50.4 

15.1 

0.00278 

0.0254 

9.1 

Bottom 

1 

50.2 

15.1 

0.00254 

0.00740 

2.9 

Middle 

1 

50.2 

15.1 

0.00238 

0.00331 

1.4 

Top 

2/3 

48.8 

14.5 

0.00171 

0.00240 

1.4 

66/19.0 

0.05/0.95 

Whole 

1 

65.8 

18.9 

0.00520 

0.0240 

4.6 

Bottom 

2/3 

48.7 

13.3 

0.00126 

0.00383 

3.0 

Middle 

2/3 

57.6 

16.2 

0.00158 

0.00189 

1.2 

Top 

1 

66.2 

19.1 

0.00216 

0.00323 

1.5 

5.  DISCUSSION 


The  aim  In  writing  TGA2  was  to  provide  the  means  of  calculating 
kinetic  parameters  from  dynamic  thermogravimetric  data,  and  to  do  it 
easily  and  quickly,  and  for  the  results  to  be  at  least  as  reliable  as 
other  available  methods. 

Doyle  [6J  has  provided  adequate  thermogravimetric  data  on 
octamethylcyclotetrasiloxane  (OMCl'S)  and  polytetrafiuoroethylene  (PTFE) . 
The  relevant  numbers  (19  temperatures  and  1 heating  rate  for  each  run) 
were  abstracted  from  the  paper  and  applied  to  the  computer  with  program 
TGA2.  The  process  for  both  sets  took  approximately  30  min  plus  a few 
seconds  computer  time.  A comparison  of  the  results  is  shown  in  table  8. 


* 


TABLE  8 


A Comparison  of  Results  from  Reference  [6]  and  TGA2 


Method 

OMCTS 

PTFE 

Activation 

Energy, 

kcal/mole 

Order 

Activation 
Energy , 
kcal/mole 

Order 

Ref  [6] 

11.65 

0 

66-68 

1 

TGA2 

11.44 

0 

60.45 

2/3 

Zsako  has  published  data  [7]  from  the  thermal  decomposition  of 
(I)  Co  (dime thy lglyoxirae)  2 (p.  ethylaniline)  2NCS  and  (II)  Co(dimethylglyoxime)2 
(y-picoline)2NCS.  He  has  calculated  kinetic  parameters  for  the  decomposition 
of  these  compounds  using  the  Horowitz  and  Metzger  method,  the  Freeman  and 
Carroll  method,  as  well  as  the  modified  Doyle  method.  These  results  are 
compared  with  TGA2  in  table  9. 


TABLE  9 


A Comparison  of  Kinetic  Parameters  Calculated  for 
Cobalt  Complexes  I and  II  by  Different  Methods 


Compound  I 

! 

Compound  II 

Method 

1 

E 

1 

Reaction 

Order 

— 

Z 

S* 

E 

1 

Reaction 

Order 

Z 

s* 

iHorowitz  and  Metzger 

42 

2 

- 

- 

33 

1 

- 

+ 1 

Zsako 

28.9 

1 

11.71 

- 5.84 

28.3 

1 

11.45 

- 7.06 

Freeman  and  Carroll 

29.2 

1.16 

- 

- 

28.7 

0.82 

- 

- 

TGA2 

27.7 

1 

10.77 

- 10.16 

28.5 

1 

11.12 

- 8.54 

Both  tables  8 and  9 indicate  that  the  TGA2  results  are  comparable 
with  the  original  results  of  Doyle,  Zsako  and  with  those  of  Freeman  and 
Carroll*  All  of  these  give  a lower  value  for  activation  energy  than  the 
Horowitz  and  Metzger  method. 


A word  of  warning  is  probably  necessary  at  this  stage. 

The  method  of  using  TGA2  is  so  simple  that  there  is  a 
temptation  to  do  the  calculations  and  quote  values  for  kinetic  parameters 
which  are  meaningless,  and  it  is  especially  important  to  remember  the 
limitations  of  the  method;  the  important  ones  are  listed  as  follows:- 
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(a)  The  weight  loss  step  which  is  considered  must  be  directly 
related  to  a simple  reaction  and  not  compounded  from  several 
different  processes. 


* 


(b)  The  results  will  only  be  as  good  as  the  original  TGA 
experimental  results  allow. 

For  instance,  care  must  be  taken  to  minimise  diffusion 
effects  which  arise  from  the  use  of  samples  with  a large 
particle  size,  or  a high  heating  rate.  Particular  attention 
should  be  paid  to  the  lower  1/3  or  "tail"  of  the  curve  which  has 
been  shown  to  be  distinctive.  Also  in  this  area,  a build-up  of 
non-volatile  catalysts  as  the  weight  loss  proceeds  may  cause 
spurious  results,  eg,  impure  carbon  burning  in  air. 

(c)  The  third  important  consideration  is  to  decide  which  mechanism 
best  describes  the  solid  state  reaction.  Zsako  has  used  the  six 
relatively  simple  equations  listed  in  table  1,  but  many  others 
have  been  developed  and  they  are  designed  to  accommodate  different 
rate  controlling  processes.  Examples  are  given  in  table  10. 

TABLE  10 


Kinetic  Eauations  Based  on  Various 

Solid  State  Mechanisms  [8] 

Function 

Equation 

Rate-Controlling  Process 

a2  “ kt 

One-dimensional  diffusion 

d2 

(1  - a)  In  (1  - a)  + a - kt 

Two-dimensional  diffusion, 
cylindrical  symmetry 

d3 

[1  - 3/l  - a]2  » kt 

Three-dimensional  diffusion, 
spherical  symmetry;  Jander 
equation 

d4 

(1  - jx)  ~ V(1  - a)2  ' kt 

Three-dimensional  diffusion, 
spherical  symmetry; 
Glnstllng-Brounshtein  equation 

Fi 

- In  (1  - a)  - kt 

Random  nucleation,  one  .nucleus 
on  each  particle 

A2 

/-  In  (1  - a)  ■ kt 

Random  nucleation;  Avrami 
equation 

a3 

V-  In  (1  - a)  - kt 

Random  nucleation;  Avrami 
equation 

r2 

1 - A - a - kt 

Phase  boundary  reaction, 
cylindrical  symmetry 

Ra 

1 - Vl  - a « kt 

Phase  boundary  reaction, 
spherical  symmetry 

Notable  methods  of  calculating  kinetic  parameters  from  thermal 
analysis  curves  have  been  developed  by  Freeman  and  Carroll  [12],  Coates 
and  Redfera  [13],  Ozawa  [14]  and  Sharp  and  Wentworth  [IS].  These  have 
been  compared  by  Ozawa  [16]. 
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APPENDIX  A 


PROGRAM  TGA2 

IMPLICIT  RE  AL  *8 ! A— H ,Q— Z > 

DIMENSION  bl 15,50) 

DIMENSION  OTDH! bO ) 

DIMENSION  EAKT(SO) 

0 1 MENTION  6115,50) 

DIMENSION  HI5Q 1 
DIMENSION  ORDER!  lb) 

DIMENSION  PX150) 

DIMENSION  T 150 ) 

DIMENSION  TClbO) 

PlANCK*6.629E-2  7 
BOLTZ*1.3eOE-lfc 
R = i .98700 
KE aD( b 1 10)  NOATA 

10  F3KMAT1I2) 

RcAO I b , 1? ) HD1F 

15  FORMAT (F6.3 ) 

KfcA0l5,.6>  (OKJEfU  J),  J=i  ,15) 

16  FORMAT UbAA) 

READ! 5,20)  !H(I ),I*i, NOATA) 

20  FOkMAT !F7.A  I 

100  RE AO (5, 8 H RUN, STEP, RANGE , QM 
3 FORMAT  1 18,2A8,F8. 3) 

0=0M/60. 

IF!IRUN).0:,10*,102 

101  STOP 

102  REAOlb ,21)  (TCII),IS1 ,NOATA) 

21  F Ok  MA  T ! F 7 • 1 > 

DO  22  1=1, NOATA 
T! 1)=TC! I I ♦273. 200 

22  CONTINUE 

CALL  SLOPE  1 T ,N0ATA,H01F,l»T0H) 

9H3T=l./0TDHil) 

WHITE  16,9) IHUN, STEP, RANGE, NOATA,QM 

9 FORMAT!*!*,  *RUU*, 18, 2A8, AX, »NDATA=*,l2, AX,  ‘HEATING  RATE  =*,F8.3,*D 
1EG.C.PER  MIN. • / / ) 

L=!N0ATA+l>/2 
WHITE  16, 17)  TLIl) 

17  FORMAT  ( * TEMPERATURE  FOR  HALF  WEIGHT  LOSS  IN  THIS  STEP  = *,F7.1,2X 
1,*0EG.C.  '/) 

GAPRQX=l  -R*|TIA)**2)*D>«T)/!-DL0GtHlA) ) | 

WHITE! 6, 13 ) H(4|,0H0T 

18  FORMAT!*  THE  SLOPE  DH/DT < AT  H=*,F7.A,*)  =*,F12.2/) 

WRITE16, *91  EAPROX 

19  FORMAT!*  A HOUGH  ESTIMATE  OF  ACT. ENERGY  =• , -3PF10.2 ,*  KCAL/MOLE* 
1/) 

WH I TE ( 6 , 11) 

11  FORMAT!*  THE  bCST  FIT  IS  INDICATED  bY  THE  MINIMUM  SIGMA  IN*/*  TH 
IE  FQLOWING  TAbLt  SO  USE  THE  VALUES  IN  THAT  LINE.*//*  NOTE, THEY  A A 
IE  EMPIRICAL  A(4U  OEPF.NU  ON  THE  E XPTL  .MF  TrlOO  . • / / } 

WHITE ! 6, 5D ) 

50  FORMAT! 9 9X, *Z*/*  REACTION* ,IbX, *ACT. ENERGY* ,12X, *PRE  EXP. *, AX ,• AC T 
1.  ENTROPY*  /3  X,  * ORDER  *,  7X  , • SI  GMA*  , 3X  , • IKC  Al/MOLE  )*,AX,*LOG 

1 Z ■ , AX  , * FaC  TOR  * ,3X,  •ICAL/DEG.NOLE)*/) 
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THIS  PAGE  IS  BEST  QUALITY  PRACtICABHi 

FROM  OOrY  TO  DDC  _ 


1 


00  56  J=  1 * 6 
ITEST-2 

5 CALL  VDOIAOI ITEST,EAPRQX,F, 100,25. 000,0.000100, 500. ODO) 

00  TO  (1  ,2,  3,4)  ,1TEST 
1 00  *4  I=I,NDATA 

110  Gll,l)=l.-H(I) 

111  G42,I  >=3./2.*tl.“(Hm**l2./3.||) 

112  G(3,l)=2.4(i. -OSQRT (H (1 ) ) ) 

113  GI4,1)»3.*I  i.-IHIl)**  (A./3.IM 

114  G(5,I  )*-0L0G4Hll  )) 

115  g!6,I)=<  1.-HIUI/HIII 

1 46  EART ( 1 ) =t APROX/ (R*T (1)1 

IF! EART! 1 ) ■"2.  IDO  >109 ,109, 117 
109  F=1000.000 
GO  TO  5 

117  PX 1 1 )*OEXP  C -E  ARK  I ) ) * 1 1 1./  I EART ( I >**2  > >— ( 2. / ( EART ( I )**3  ) ) ) 

118  B<  J,  I > = l OLOGIGt  J,1>)-0L0G!PXI  im/2.30300 
4«,  CONTINUE 

F=0.000 

FS=0.000 

DO  45  1=1 ,NDATA 

119  FsF*b(J,  1>4*2 

120  FS=FS-*B!J,I) 

45  CONTINUE 


F=QSQRT  ( !F- 1 FS**2  )/FL OAK  NOAT A ) )/  FLOAT  INOAT  A-I ) ) 

GO  TO  5 

2 I F| EART INDATA)— 4. 000)252,251,251 

251  BBAR=FS/FLOAT (NOATA ) 

2LOG=BBAR»OLOGIR40)/2.303DO-OLOG(EAPROX|/2. 30300 
Z=10.**2L0G 

S=R*DL0GKZ*PLANCK)/4B0LT2*TIL))> 

WRITE  I 6, 60)  ORDER! J),F,EAPKUX,2LUG,Z,S 
60  FORMAT  (JX,A4,3X  , E12 .3  ,— 3PF10.2 , 0PF10.  2,E12.3,F13.2) 

GO  TO  55 

25 2 WRITE  46,253)  ORDER! J) ,F,EAPROX 

253  FURMAT !3X,A4,3X,E12 .3 ,“3PF 10. 2, • CALCULATION  INACCURATE  AT  THIS  LO 
1W  VALUE.  ') 

GO  TO  55 

3 WRITE  I 6, 65 ) 

65  FORMA T ( 20X , 'ROUNDING  ERRORS*) 

GO  TO  55 

4 WRITE (6,70) 

70  FORMAT (20X, *MAXFUN  VALUE  TOO  SMALL*//) 

55  CONTINUE 
GO  TO  100 
ENO 


^»s“‘ 


s 


1 


c 


c 

c 


SUBROUTINE  SLOPE (Y,N,XD1F,DYDX> 
IMPLICIT  KEAL*d(*-H,0-Z) 
DIMENSION  0V0XI50) 

DIMENSION  Y ( 50 1 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


□YF I C 50) 
DYF/( 50) 
DYF3I50) 
DYFM  50  ) 
DYF5( 50) 
YMUii 50) 
YMU3I50) 
YMU5I50) 


NMl=N-i 
DO  5 1=1, NMi 
OYPlll)=Yl 
5 CONTINUE 
NM2=N-2 
DU  10  1=1, NM2 

YMUl(I)=CUYFUl)+0YFllI*l))/2. 

0YF2I I)=DYF1(1*I >— OYFil 1 ) 

10  CONTINUE 
NM3=N-3 
DO  1.5  l = i.,NM3 
DYF31I )=DYF2< 1*1 I-0YF2I 1) 

15  CONTINUE 
NM5=N-5 
DO  20  1=1, NM5 

YMU3 ( 1)=(DYF3(1 ) +DYF3 tl+l))/2* 

DYFM  1)=DYF31I*1  I-UYF31 1) 

20  CONTINUE 
NM5=N— 5 
DU  30  1=1, NM5 
DYF5<  I)=DYF4<I*1)-DYFM1> 

30  CONTINUE 
NMo=M-6 
DU  35  I*i,MM6 

YMU5U)=IOYF5ll)-DYF5tl  + i)  )/2. 

35  CONTINUE 

DO  50  1=1, NM6 

JYCX(  I ) = ( YMU1  »!♦*)-  t YMUSt  !♦!)/&.  )♦!  YM\J5(l  )/s0.  ) I/X01F 


50  CONTINUE 
RETURN 

END 


1 

Y*  >» 
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FIGURE  2 OUTPUT  FROM  PROGRAM  TGA  2 
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FIGURE  3 Effect  on  Sigma  Values  of  Progressive  Removal 
of  Top  of  Standard  Curve 
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Abstract  zsako  has  modified  a curve  fitting  method  reported  by  Doyle  and  in  this  report  the 

method  is  further  simplified  to  allow  the  use  of  a digital  computer.  The  best  description 
of  the  thermogram  step  was  found  to  be  given  by  19  equally  spaced  co-ordinates  between  the 
(1  - a)  limits  of  0.05  to  0.95.  These  co-ordinates,  together  with  the  heating  rate,  are 
used  to  calculate^ E and  log  Z,  by  an  iterative  process  which  minimises  the  standard  devia- 
tion of  the  log  g(a)  - log  p(x)  values.  The  standard  deviation  allows  a comparative  assess- 
ment of  the  validity  of  the  kinetic  equations  chosen  to  represent  the  mechanism  of  the 

reaction.  A further  simplification  involved  the  use  of  an  approximate  equation  for  - log 

p(x)  in  place  of  tabulated  values.  The  TGA2  method  of  calculation  was  shown  to  be  valid  by 

re-calculating  E and  log  Z for  standard  curves  which  were  generated  from  tables  of  - 
log  p(x).  Published  data  by  Doyle  and  Zsako  were  re-calculated  and  the  answers  were  in 
good  agreement  with  the  original  results.  
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So—  Metric  and  SI  Unit  Conversion  Factors 


(Based  on  DEF  STAN  00-11/2  "Metric  Unit*  for  Ua*  by  the  Ministry  of  Defence", 
DS  Met  5501  "AWRF.  Metric  Guide"  and  other  British  Standards) 


Quantity 

Unit 

Syabol 

Conversion 

Basic  Units 

Length 

■atrt 

m 

In-  3.2808  ft 

1 ft  - 0.3048  n 

Mass 

kllogr— 

1 kg  - 2.2046  lb 

1 lb  - 0.45359237  kg 

1 ton  * 1016.05  kg 

Derived  Units 

Force 

nevton 

N - kg  n/s2 

IN*  0.2248  lbf 

1 lbf  - 4.44822  N 

Work,  Energy,  Quantity  of  Heat 

joule 

J • N 

1 J - 0.737562  ft  lbf 

1 J - 9.47817  x 10~4  Btu 

1 J - 2.38846  x 10~4  kcal 

1 ft  lbf  - 1.35582  J 

1 Btu  - 1055.06  J 

1 kcal  - 4186.8  J 

Power 

vatt 

W - J/s 

1 V - 0.238846  cal/s 

1 cal/s  - 4.1868  W 

Electric  Charge 

couloab 

C - A a 

- 

Electric  Potential 

volt 

V - V/A  » J/C 

- 

Elactrical  Capacitance 

farad 

F - A s/V  - C/V 

Electric  Resistance 

ohm 

n - V/A 

s • i a'1 

Conductance 

si— en 

Magnetic  Flux 

weber 

Wb  - V s 

Magnetic  Flux  Density 

tesla 

T - Vb/n2 

- 

Inductance 

henry 

H - V s/A  • Vb/A 

“ 

Co-lax  Derived  Units 

Angular  Valoclty 

radian  per  second 

rad/s 

1 rad/s  - 0.159155  rev/s 

1 rev/s  » 6.28319  rad/s 

Acceleration 

netre  per  square  second 

n/s2 

1 n/s2  - 3.28084  ft/s2 

1 ft/s2  - 0.3048  n/s2 

Angular  Acceleration 

radian  per  square  second 

rad/s2 

- 

Pressure 

newton  per  square  — tre 

N/n2  - Pa 

1 N/n2  - 145.038  x 10-6  lbf/in? 

1 lbf /in?  - 6.89476  * 105  N/n2 

bar 

bar  - 10s  N/n2 

1 In.  Hg  - 3386.39  N/n2 

Torque 

nevton  netre 

N a 

1 N n - 0.737562  lbf  ft 

1 lbf  ft  - 1.35582  N m 

Surface  Tension 

newton  per  —tre 

N/n 

1 K/m  - 0.0685  Ibf/ft 

1 lbf/ft  - 14.5939  N/n 

Dynaalc  Viscosity 

newton  second  per  square  —tre 

N a/a2 

1 N s/n2  - 0.0208854  lbf  s/ft2 

1 lbf  s/ft2  - 47.8803  N s/n2 

Kin— atic  viscosity 

square  netre  per  second 

n2/a 

1 n2/s  - 10.7630  f t2 / s 

1 f t2/s  - 0.0929  n2/s 

Ther— 1 Conductivity 

watt  per  —tre  kelvln 

W/n  K 

- 

Odd  Units* 

Radioactivity 

becquerel 

Bq 

1 Bq  - 2.7027  x 10_n  Cl 

1 Cl  • 3.700  x 1010  Bq 

Absorbed  Dose 

gr»y 

Gy 

1 Gy  ■ 100  rad 

1 rad  - 0.01  Gy 

Dose  Equivalent 

s lever t 

Sv 

1 Sv  - 100  r— 

1 r—  - 0.01  Sv 

Exposure 

couloab  per  kllogr— 

C/kg 

1 C/kg  - 3876  R 

1 R - 2.58  » 10-4  C/kg 

Rate  of  Leak  (Vacuua  Syst— s) 

allllbar  litre  per  second 

nb  1/s 

1 nb  - 0.750062  torr 

1 torr  - 1.33322  nb 

♦The**  tens  are  recognised  tens  within  the  —trie  syst— . 


